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The characteristics of the star model with the "prescribed" density distribution p — p c [l — (r/R) a ] 
are analytically studied. The model has been first considered briefly in our 30-year old note of a 
restricted circulation Jl5[ . 

Introduction 

We choose the distribution of density p(r) in a spherical-symmetric star in the form 

p(r)=p c (l-x a ), x = (1) 

here p c — p(0) is the central density, r is a running radius, R is the radius of the model, a is a free parameter, (a > 0, 
as we require dp/dr < 0). 

Central-to-mean density ratio 
From Eq. (|l|) we have for a cumulative mass m, total mass M and for the central-to-mean density ratio p c /~P'- 

m ( x )=A7rp c R 3 (^--^); M = m(l) = ±* Pc R 3 ^ - 4f pi? 3 ; f = *±*. (2) 



Central moment of inertia 

The moment of inertia of the star about its center is: 

1= r*d m =^p c R 5 -^: (3) 
Jo 5 5 + a 

with help of Eq. (2) we may write down I in terms of M and R: 
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Gravitational potential energy 

The total potential energy of the star is: 

w _ n [ M m J m _ 16tt 2 a 2 (ll+2a) p r>5 n 2 _ -3 (3+a) (11+2 a) G M 2 (<?\ 
VV — U JQ r arn - 15 (3+a)(5+ii)(5+2«) UJt P c ~ 5 (5+a) (5+2 a) R ' 

WUM-ratio 

The central gravitational potential U(0) and the potential at the surface U(R) of the spherically-symmetric star 
are 

pR r< yr 

U{0) = -4nG I p(r)rdr; U(R) = — — . (6) 
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for the model (Fll): 
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Combining Eqs (2), (5), (7) we obtain the WUM-x&tio (Of for our model: 



WUM = ^— = 2 ( 2 + Q H 11 + 2a ) (8) 
U(0)M 5 (5 + a) (5 + 2a) { > 



Pressure-density relation 



Now we find the relation between pressure P and density p in the star, which has the density distribution ([!]) at 
hydrostatic equilibrium. In other words, we find the equation of state of matter which leads to distribution of density 
(0) in the star. We remind that a star with the equation of state P(p) which does not include a temperature is called 
pseudopolytrope Q. 

The equation of hydrostatic equilibrium of spherically-symmetric star reads: 



1 dP m 

p dr r 2 



(9) 



Using Eqs (1, 2) and integrating Eq. (g) with initial condition P(0) — P c (central pressure) we obtain the distribution 
of pressure and density within the star: 



P(x)=P c -GirR 2 p c 
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3(2+q) (3+q) + (3+q) (1+q) 



p(x) = p c (l - x a ); x = r/R. 
From Eq. (0) using condition P(l) = 0, we obtain the following relations for the central pressure: 



2 7ra 2 (4 + a) 
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3 (1 + a) (2 + a) (3 + a) rc 8 tt (1 + a) (2 + a) R 4 ' 
From Eqs. (0, [□]) we obtain the following relation between P and p ("equation of state", EOS): 



£- = l 

Pc 



+ 1 1 _ yj i + a y + yy ) 



2 p_ i 3(2+q) /_p_N 2 



(10) 
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Note that for any given pair of parameters p c , P c , (and fixed a) we have a particular EOS given by Eq. (12). 
Therefore these particular pseudopolytropes have the two-parametric EOS, while e.g. the classic polytropes with given 
adiabatic index 7 have one-parametric EOS: P = K p 1 ', the parameter K being related with polytropic temperature, 
see for more details If we calculate a model with other values of central density p* and pressure P* with the 

same EOS (12), the resulted distribution of density will not coincide with the law (|], see the next Section). Particular 
cases a = 1 and a = 2 give the more simple equations of state: 
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(13) 



f = ^) 2 (i + ^); a = 2. 

-Tc 2 p c p c 



(14) 
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FIG. 1. Equations of state for different values of a: 



a = 1, 



with 7 = 4/3. Abscissae are values of density ordinates are values of pressure P/P c 



a — 1/2, and solid line - polytropic EOS 
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Although these equations of state, and the more general EOS @ are derived and are formally valid only for 
— < 1, we will loosely use them also at larger densities, see the next Section. 



Density distribution in the a = 2 model 

Here we calculate density distribution of the particular model with a = 2, which corresponds to EOS given by 
Eq. (|l4|). With this EOS, the equation of hydrostatic equilibrium (||) is reduced to the next equation for the density 
distribution p(x): 

sf<!> + i>3Sl"««V (15) 

In deriving this equation, we used the first relation in Eq. (11). Important notice is that Eq. (|l5|), at the case 
p(0) — 1, has the solution which corresponds to the law (Q) (with a = 2), and all parameters corresponding to this 
case are taken to be equal to unit (that is, at p(Q) — 1, we have R = 1, M = 1 and so on. Also, the dimensionless 
running radius in Eq. ( |l5[ ) x is expressed in units of the total radius of star R, which is not equal to 1 at p(0) ^ 1. 




FIG. 2. Density distribution for EOS P = p ^ +1 ' > a ^ three values of central density p(0) = 2, 1 and 1/2. Abscissae are 

pO) 
p(0) ■ 



values of dimensionless running radius x = r/R, ordinates are values of 



Eq. ( |15|) was solved numerically for two values of p(0), 2 and 1/2. In these two cases the radius of the star is 1.169156 
and .9070496, respectively. The resulting density distribution is shown in Fig. 2, where abscissae are x — r/R and 
ordinates are p(x)/p(0); p(0) = 2 case corresponds to the upper solid curve, while p(0) = 1/2 corresponds to the 
lower solid line. Also shown is the standard distribution - Eq. (Q) with a — 2, dash line in Fig. 2. For values of the 
mean-to-central density ratio 7j/p(0) we have: .363325, .4 and .442107 at p(0) equal to 1/2, 1 and 2, respectively. For 
larger central density, EOS (|14|) has larger adiabatic index 7, EOS is more stiff, and this leads to more homogeneous 
configuration with larger mean-to-central density ratio. 

Pade approximants 

We briefly describe the Pade approximants method which we use further in this paper. 
From the segment of the series: 

fe=4 

s i = l + Y j a k x k , (16) 
fc=i 

coefficients of Pade(2,2) approximants H 

1 + A x + B x 2 , , 

PA = l + Cx + Dxi ^ 

are defined as follows Q: 

a| - a 2 q 4 ai a 4 - Q2 Q3 n . n , n A , n nsA 

V = —5 r5 C = — 5 B — a2 + ai G + D\ A = a\ + G. (18) 

a^ — a\ a A a^ — a\ a 6 
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Pade approximants to density distribution 



From Eq. @ we have series expansion at the center: 

- 5 . . _ -15(-1+P(0)) . 

2+3p(0V " 2 C2+3of0)) 3 ' 

(19) 



_ 150(l-9p(0)) (-l+p(0)) , _ 25(-l+p(0)) (-5+222p(0)-837p(0) 2 ) 

° 3 7(2+3p(0)) B ' 04 — 7(2+3p(0)) 7 

Evidently, at p(0) = 1 the series is reduced to p(x) = 1 — x 2 . Expressions for coefficients of Pade(2,2) approximants 
are too cumbersome at the general case of arbitrary p(0) so we present only one particular case: 

„ , . „ 12816384- 13840960 x 2 + 3314065 x 4 , nn . 

p(0)=2: PPade {x)=2 12816384 _ 5830720^+45345^4 • ( 20 ) 

Using this "analytic solution" we find the radius of the model equal to 1.17720 which is 1.0068 times the calculated 
value 1.169156. Also, we calculated mean-to-central density ratio according to ( p0| ) and found value 0.435602 which 
is 0.98528681 of the numerical value .442107. 



The adiabatic index is defined as: 



Mean adiabatic index 



d In P 



d In p 



Distribution of adiabatic index within the model is: 



7(x) 



2 (1 + a) (2 + a) x 2 ~ a (3 + a - 3 x a ) (1 - x a Y 



[4 a 2 + a 3 - (6 + 11 a + 6 a 2 + a 3 ) x 2 + (12 + 14 a + 2 a 2 ) x 2+a - (6 + 3 a) x 2 + 2 a ] ' 



(21) 



(22) 



At the surface of the star (at x = 1) 7 = 2 for any a while the value of 7 at the center of the star depends on a: 
7(0) = at a < 2, 7(0) = 00 at a > 2, and 7(0) = 5/2 at a = 2. 

Moreover the function 7(2;) is non-monotonic at a < 2, while at a > 2 the function 7(2) is monotonically decreasing 
with increasing x, relative running radius, see Fig. 3. 
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FIG. 3. Adiabatic index as function of x — r/R, according to Eq. (22), curves from the lowest one to the upper one 
correspond to values of a = 1/2, 1, 2, 2.1, respectively. 

Therefore at a < 2, in the inner parts of the star, the value of 7 is less than 4/3, the critical value of adiabatic 
index, while in the outer regions 7 is always larger than 4/3. 

The stability of the star against radial perturbations is defined by the mean value of adiabatic index. That is why it 
is of interest to calculate the mean adiabatic index, 7, which is defined as follows: 



1 



7 P dv; 



Iv Pdv Jv 

here an integration is over the volume of star. We obtain at a < 5: 

12(7 + a-) 
(5-a)(ll + 2a)' 



7 = 



(23) 



(24) 



For any a > 5, 7 — > 00 due to divergence of 7(7-) P(r) r 2 at the center of the star. The function 7(a) is monotonically 
increasing and even at a — 0, 7 = 84/55 > 4/3, that is the model is always stable against to the radial perturbations. 
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Ellipticity of slow rotating star 



For the slow rotation, the ellipticity distribution inside the spherically-symmetric star is governed by the Clairaut's 
equation (|]: 

[ e^)_6e^) ] f X pit)t 2 dt + 2p{x)[ eM +e > ix)]=0 . (25) 

X X J q X 

Here x — r/R is dimensionless running radius, < x < 1, and e(x) is an ellipticity of the equidensity surfaces 
within the star, e(x) = 1 — r p (x)/r eq , with r p and r eq being the polar and equatorial radii of the equidensity surfaces 
which are assumed to be the biaxial ellipsoids of revolution. 

The Equation (E5ft can be solved in terms of the series expansion at the center, which we write down here for two 
cases: 



a = 


1 




e(x) 


— 2^i=0 


Cj X , 


_ _ 3 i 2 +5i-4 . 
L * 4 i(i+5) 


i > 1; 


e(x) 
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e(x) 




Ci X , 


3 2i 2 +9i+2 
t-i+1 — 5 (2i+7) 


Cj_i i > 0; 


e(x) 
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,6 a; 2 
r 35 


, 13 x 4 
"r 175 


i 52 x 6 
' 1375 


, 141 x 8 , 1974 x 10 
' 6875 1 171875 ' 





(26) 



We solved Eq. (]25|) numerically for the model (1) with two values of a, see Figs. 4 and 5. 
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FIG. 4. Distribution of ellipticity within the star for a = 1. Solid line - exact values of e(x), dash line - series expansion 
according to Eq. (26). Abscissae are values of running relative equatorial radius x — r eq /R eq , ordinates are values of ellipticity 
e(x). 
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FIG. 5. Distribution of ellipticity within the star for a = 2. See caption to Fig. 4. 

The values of the surface-to-central ellipticity ratio of ellipticity e(l)/e(0) are 1.6608 and 2.0488 at a equal to 1 and 
2, respectively. In the theory of (slowly) rotating stars the ratio of centrifugal acceleration w 2 to the gravitational 

acceleration GM/R at the equator, m = ^rjf is introduced. Also, the ratio m-to-e(R) ratio is introduced, which 

in the terms of functions e(x) and e'(x) is ^ = | (2 + ^jy)- F° r the model in question, we have a — 1 : e(l) = 

1.6608, e'(l) = 2.0488, m/e = 1.2936, and a = 2 : e(l) = 1.3312, e'(l) = 1.3785, m/e = 1.2142. 



Pade Approximants for the ellipticity distribution 



For the ellipticity distribution inside the star, using the series expansion (26), we have Pade(2,2) approximant 

, s 5 (25025 - 14036 x+ 1138 x 2 ) , s 

a = 2 : e Pade (x) = —7 rf. (27) 

padev ; 7 (17875 - 13090 x+ 1729a; 2 ) v ' 
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Even at the surface this expression is very accurate: ep at 
solution for the Clairaut's equation for a = 2 model. 
For the a — 1 model Pade(2,2) is also rather accurate: 



J (l)/e ca ; c (l) = 0.99891. In fact (§7|) gives the analytical 



672 - 448 x + 41 x 2 
7 (96 - 88x+ 15x 2 ) 



(28) 



With epade(l)/e ca lc(l) = 0.99108. 



Rotation, contraction and mass loss 



The abovementioned analytical expressions for parameters of star may be applied to study of evolution of rotating 
star, or to evaluation of pulsational periods of star or to any other problem where the simplicity and explicity of 
evaluation justify the approximation certain inaccuracy. 

We consider in this section the evolution of contracting and rotating star with mass loss. Let the rotating star be 
contracting in such a way that the distribution of density follows the law ([!]). At the some moment the condition 
for mass loss may be reached at the star's equator. We assume that the star is spherically-symmetric and rotation is 
steady-state (that is angular rotation is constant over the volume of star, uj{t) — constant) (that may occur at slow 
contraction and rapid rotation), then we may write the mass loss condition as: 

? G M 

" - -w- (29) 

The rate of angular momentum loss is dL — ujR 2 dM, and taking into account that L = kwM R 2 , (k is a structure- 
dependent parameter) we have: 

M f 



M -\R J • ^—k- ^ 

As a result, a value of the final mass of the star depends only on parameter k, (see also The values of fc for 
polytropes, white dwarfs and "stepenars" are given in [5 - 8]. For the star with the density distribution ([!]): 

5 (5 + a) 



note that in Eq. (4) the moment is central, and axial moment is I axis — 1/3 1 center- 
At a = .06 that corresponds, by value of ratio — , to n — 3 polytrope, we have k ~ .012 and /3 ~ .078. Rotation- 
induced distortion of figure from sphere leads to (3 ~ .02 pjj. For the star in the pre- white dwarf stage, a value of (3 
should be even less as the star in pre-white dwarf or pre-neutron star stage is a (very) hot star with elusive extent 
envelope and with small dense core (and with very small value of k). The envelope contains the large rotational 
momentum and the small mass therefore the mass loss per unit momentum loss is very small for a such star. 
At the other hand for the homogeneous star, k = .2, /3 = 1/7 that is 7 times the value of n — 3 polytrope. Therefore 
if the pre-white dwarf star would be more homogeneous then at contraction the rotating star could lose more mass. 
However even in this case to reach a sizable value of mass loss at the pre-white dwarf stage, the star should have the 
rotational momentum by order of magnitude larger that the main sequence star, and for neutron star the difficulty is 
even larger. The presence of factors leading to the more homogeneous structure could lead to reducing the difficulties 
in explaining the origin of white dwarfs and neutron stars by the mass loss from ordinary stars. 



Pulsational periods of pre-white dwarf stars 



Ledoux and Pekeris EJ] using the energy method obtained the following formula for the frequency of the fundamental 
mode of the adiabatic radial pulsations of the spherically-symmetric star: 



-W 

<r 2 = (3 7 - 4) — . (32) 



Using the formulas (4, 5, 12) we obtain 



., 32 7T (l + a)(4 + a) 



3 (5 — a) (5 + 2 a) 
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At p ~ 10 4 (a star with solar mass and radius ~ 1/20 the solar radius) and a = 1, we have a ~ 8.9 10 2 , or for the 
pulsational period P — 2 n/a ~ 70s. 

Such and even larger periods of light variations occur in HL Tau stars, G 44-32 and R 548 [jl(| , (ll) which are apparently 
at the late s tag es of their evolution, and probably - at pre-white dwarf stage. Recent discussion of pre-white dwarf 
stars see in 14|. 



Radial pulsations 

The differential equation for the adiabatic small radial pulsations of the spherically-symmetric star is: 

t 7 P S) + y 9 ^ + Tr [(3 7 " 4) P] } = ° ; (34) 

here y = Sr/r is the ratio of the radial displacement to the radius, and a is the frequency of the pulsations. At 
the center (at r = 0), we have condition 6r — ry — 0, at the surface, r — R, we require that amplitude is finite. We 
consider here the case of a — 2 when Eq. ( |34| ) is reduced to the form: 

2 y{x) x (19 - E 2 - 15 x 2 ) + 2 (-10 + 29 x 2 - 15 x 4 ) + y"{x) (-5 x + 8 x 3 - 3 x 5 ) = 0; (35) 

here x — r/R, and E 2 is dimensionless value, expressed in units of y-jjj — ^j- G~p, see Eqs. (2,11). 



Equation (35) can be solved in terms of series expansion, and the recurrence relation can be written down for 
coefficients of the power series which may contain only even powers of x. Note that Eq. ( |35| ) is the linear differential 
equation and the function y is defined up to an arbitrary factor, so we may take y(0) = 1. Then the first terms of the 
series expansion at the center are: 

19 - E 2 , 490 - 104E 2 + E 4 , 2670 - 12610 E 2 + 239 E 4 - E 6 fi , . 

y = l-\ x 2 H x 4 H x 6 . 36) 

y 25 1000 60000 v ' 

At the surface (x = 1) we have the series expansion 

»r ,E , ll2 2E ~\~ E , . On 

y(x) = y(l) [1 + (— - 1) (1 - x) + (1 - x) 2 ]. 

For arbitrary x the Eq. ( |35| ) may be solved numerically. We calculated the four lower modes of radial pulsations 
and found for the corresponding eigenvalues E 2 the values 10.325, 39.083, 81.04, and 136.1, see Fig. 6. 




FIG. 6. Eigenfunctions for 4 lower modes of radial adiabatic pulsations as functions of x = r/R for the star model with 
a = 2, see Eq. (35). Curves from the lowest (fundamental) mode (no nodes) to upper mode (3 nodes) correspond to values of 
E 2 equal to 10.325, 39.083, 81.04 and 136.1, respectively . 

The lowest eigenvalue corresponds to a 2 = 10.325 * 2/3ttG7j = 6.88 7tG7j =, that only slightly differs from value 
obtained by energy method, 647t/9G7j = 7.l7rGp, see Eq. ( p3[ ) with a = 2. Note that the evaluation of a given by 
Eqs. (B3, p3h for the frequency of the fundamental mode is to be either equal to or greater than the true value M . 



The evaluation of pulsational frequency of the lowest mode can also be obtained directly from Eq. (35), if we put 



y = const (which is not too rough assumption for the fundamental mode eigenfunction, see Fig. 6), then we have 

E 2 = 19-15x 2 . (37) 
We can take averaged value of right side of this equation with weight x 2 : 



1 r 1 



x 2 dx Jo 

which is very close to the "exact value" 10.325. 



(19- 15x 2 )x 2 rfx = 10, (38) 
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Conclusion 



In conclusion, we present here the analytical model density distribution which allows to evaluate qualitatively and 
quantitavely the many important characteristics of the star including the ellipticity distribution within the rotating 
star, the pulsational periods and the mass loss at the contraction stages of evolution. The certain inaccuracy of 
the approximation is the modest price paid for the large simplicity and explicity of the model. Surely, as the first 
approximation or at least as the pedagogical tool the model is of some interest. 
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